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ABSTRACT 

Star formation and accretion onto supermassive black holes in the nuclei of galaxies are the two 
most energetic processes in the Universe, producing the bulk of the observed emission throughout its 
history. We simulated the luminosity functions of star-forming and active galaxies for spectral lines 
that are thought to be good spectroscopic tracers of either phenomenon, as a function of redshift. We 
focused on the infrared (IR) and sub-millimeter domains, where the effects of dust obscuration are 
minimal. Using three different and independent theoretical models for galaxy formation and evolution, 
constrained by multi-wavelength luminosity functions, we computed the number of star-forming and 
active galaxies per IR luminosity and redshift bin. We converted the continuum luminosity counts into 
spectral line counts using relationships that we calibrated on mid- and far-IR spectroscopic surveys 
of galaxies in the local universe. Our results demonstrate that future facilities optimized for survey- 
mode observations, i.e., the Space Infrared Telescope for Cosmology and Astrophysics (SPICA) and 
the Cerro Chajnantor Atacama Telescope (CCAT), will be able to observe thousands of z>l galaxies 
in key fine-structure lines, e.g., [Si II], [O I], [O III], [C II], in a half-square-degree survey, with one 
hour integration time per field of view. Fainter lines such as [O IV], [Ne V] and II2 (O-O)Sl will be 
observed in several tens of bright galaxies at 1 < z < 2, while diagnostic diagrams of active- nucleus vs 
star- formation activity will be feasible even for normal z^l galaxies. We discuss the new parameter 
space that these future telescopes will cover and that strongly motivate their construction. 

Subject headings: Galaxies: evolution, active, starburst, Seyfert - Techniques: imaging spectroscopy 



1. INTRODUCTION 

Tremendous progress in infrared (IR) astronomy was 
made in the last decade, largely driven by the launch of 
the Spitzer and Herschel space telescopes. Nonetheless, 
several basic questions remain unanswered in the fields of 
observational cosmology and galaxy formation and evo- 
lution, emphasizing the need for future missions. These 
regard our incomplete knowledge on how the primordial 
gas collapses to form new stars, on the different modes of 
star formation, and on the potentially coeval growth of 
black holes and galaxies. Below, we elaborate on each of 
these turn. We argue that fine-structure and molecular 
lines observable in IR lines can help us address them, and 
we present a set of simulations that support the need for 
future survey-oriented facilities that can make a strong 
impact on studies of galaxy evolution. 

The collapse of the primordial gas at very high redshifts 
that led to the formation of the first stars and galaxies is 
thought to have occured via H2 line emission, which acts 
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as a very effective cooling mechanism in low-metallicity 
and low- temperature (<10'* K) environments (Wise & 
Abel 2007; Obreschkow & Rawlings 2009). Even though 
II2 remains a main coolant soon after the epoch of reion- 
ization, the direct detection of H2 gas at very high z is 
yet to be achieved observational point of view. 

Recent studies of the global infrared continuum and 
molecular gas properties of galaxies in the local and in- 
termediate/high redshift Universe suggest that mergers 
and non- or weakly-interacting star-forming galaxies fol- 
low two separate Kennicutt-Schmidt relations with simi- 
lar exponents, but different normalizations (Genzel et al. 
2010; Daddi et al. 2010; Gracia-Carpio et al. 2011). 
At high-z, both galaxy mergers and accretion of cold gas 
via cooling fiows have been suggested as sufficient mecha- 
nisms to produce IR luminosities >10^^ Lq (Powell et al. 
2011), unlike in the local Universe where ultraluminous 
IR galaxies (ULIRGs) are predominantly associated with 
mergers of comparable mass spirals (Dasyra et al. 2006). 
The parameters that determine the mode of star forma- 
tion and control its efficiency are not yet understood. 
While Spitzer and Herschel helped us to identify and de- 
termine the star-formation rates of ULIRGs at z<3, the 
bulk of star formation at these redshifts is thought to take 
place in galaxies of lower luminosity (Perez-Gonzalez et 
al. 2005; Reddy et al. 2010; Ly et al. 2011). Char- 
acterizing the star formation in such systems is essential 
for understanding the formation of present-day ellipticals 
and downsizing (Cowie et al. 1996). 

Studies of local massive galaxies and theoretical mod- 
els suggest that most spheroids host massive black holes 
(Richstone et al. 1998; Fabian 1999). The observed cor- 
relations between the masses of these black holes and the 
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luminosities and stellar velocity dispersions of their host 
galaxies (Kormendy & Richstone 1992; Magorrian et al. 
1998; Ferrarese & Merrit 2000; Tremaine et al. 2002; 
Ferrarese & Ford 2000; Shankar et al. 2009; Gultekin 
et al. 2009) are remarkable, given the vastly different 
scales that they involve. The enormous difference be- 
tween the black hole Schwarzschild radius and the char- 
acteristic radius of the bulge indicates that these rela- 
tions possibly reflect the coeval formation of the two in 
a common gravitational potential. To date, this scenario 
has been primarily tested through the comparison of the 
star-formation rate vs the black hole accretion rate as a 
function of look-back time (Treu et al. 2004; Marconi et 
al. 2004; Merloni et al. 2004; Shim et al. 2009; Grup- 
pioni et al. 2011). Indeed, both activities have been 
found to peak at comparable redshifts, between 2; ~ 1- 
3. To further \mravel and understand the relationship 
between black hole growth and bulge formation we need 
to compare the shape of the mass functions of galaxies 
with those of actively accreting black holes (Dasyra et 
al. 2011). Mass estimates of black holes in obscured 
systems will be needed for this purpose. 

The complex relationship between star formation and 
active galactic nucleus (AGN) activity, which could ex- 
tend over long duty cycles (of a few hundreds of Myrs), 
constrains the quantities that can be used as pure diag- 
nostics of either activity. Present cosmological surveys 
are hampered when it comes to disentangling AGN from 
starburst activity. For example, an unambiguous AGN 
indicator in the IR is the continuum emission from dust 
just below sublimation temperature in the AGN torus 
(Schweitzer et al. 2006; Mor et al. 2009; Mullaney 
et al. 2010). However, its detection is often uncer- 
tain. The use of a torus model to reproduce the ob- 
served spectral energy distribution (SED) of galaxies can 
depend on the choice of the star-forming galaxy tem- 
plate that is simultaneously used for the same purpose. 
Selecting high-z sources with unambiguous evidence for 
hot dust-that is independent of the chosen star-forming 
galaxy template-can lead to small, biased samples of lu- 
minous AGN. Less strict criteria can lead to incomplete 
AGN samples that are contaminated by starburst galax- 
ies. This is frequently the case for AGN selected based 
on mid-IR color-color diagrams (Lacy et al. 2004; Stern 
et al. 2005; Barmby et al. 2006). The other unambigu- 
ous signature of a hard, AGN-related radiation field is 
the line emission from ionic species that require >100eV 
for their creation. A line often used for this purpose 
is [NeV] at 14.32 /im, which however was only seen in 
stacked spectra of z >1 galaxies with Spitzer (Dasyra et 
al. 2009). 

Observing new FIR-selected galaxy samples, or sam- 
pling the peak of the IR SEDs of presently known sam- 
ples, will not only enable us to address such open ques- 
tions, but to also obtain a more coherent view of high-z 
galaxy populations. The thousands of 24 /im -selected, 
IR-bright galaxies discovered by Spitzer (Papovich et al. 
2004; Rigby et al. 2004; Fadda et al. 2004; LeFloc'h 
et al. 2004; Houck et al. 2005; Huang et al. 2009), 
the hundreds of mm galaxies (e.g. Small et al. 1997; 
Hughes et al. 1998; Barger et al. 1998; Scott et al. 
2002; Borys et al. 2003; Webb et al. 2003; Coppin et 
al. 2006), (e.g. Greve et al. 2004; Laurent et al. 2005; 
Scott et al. 2008; Perera et al. 2008; Austermann et al. 



2010), the tens of thousands far-IR and sub-mm galaxies 
detected at high-redshifts by the Herschel observatory 
(e.g. Berta et al. 2010; Oliver et al. 2010), and the 
thousands of near IR dropouts (e.g. Caputi et al. 2005; 
Daddi et al. 2005; Papovich et al. 2006) known to date 
are often treated as physically disconnected populations. 

Spectroscopic surveys with future missions or facilities 
will help us address the above-mentioned points. Ex- 
pected to be at least one order of magnitude more sen- 
sitive than their predecessors, future spectrographs will 
observe ionic and molecular gas lines, in addition to the 
dust broadband features and their underlying continua, 
for large galaxy samples over wide fields of view. In the 
spirit of demonstrating what can be done in the future in 
the IR/submm domain in the field of galaxy evolution, we 
present simulations of the number of galaxies that will be 
observed in key IR lines and features at each z range with 
the SPace IR telescope for Cosmology and Astrophysics 
(SPICA) and the Cerro Chajnantor Atacama Telescope 
(CCAT). We focus on SPICA and CCAT as these are 
conceived to be large blind survey machines, unlike the 
James Webb Space Telescope (JWST) and the Atacama 
Large Millimeter Array (ALMA), that are optimized for 
deep single-source observations. 

The paper is organized as follows. In Section 2 
we present the simulations, and the physical assump- 
tions that go into them, including the adopted line-to- 
bolometric luminosity conversion functions that we mea- 
sured from observations of local galaxies. The main re- 
sults of the simulations, i.e., the number of galaxies and 
AGN that can be detected for each ionized or molecular 
gas line with z for a given flux limit, are presented in Sec- 
tion 3, followed by a discussion in Section 4. Throughout 
the paper we will adopt the standard cosmological model 
with riM^ 0.27, 17a=0.73, Ho==71 km s^^ Mpc'^. 

2. SIMULATING SPECTROSCOPIC SURVEYS IN THE IR 

The computation of the mmibcir of objects that future 
missions will detect in each linc^ was a two-fold process. 
First, we used three galaxy evolution models that are 
based on galaxy counts and luminosity functions in sev- 
eral bands of the Spitzer, AKARI and Herschel missions 
to predict the integrated 8-1000 /im IR luminosity func- 
tion at different z bins. We have considered three dif- 
ferent models with the aim of basing our conclusions on 
solid results, and at the same time quantifying the range 
of the predictions. We then computed correlations be- 
tween the IR luminosity and the luminosities of various 
fine-structure lines, molecular lines, and dust emission 
features, using samples of local galaxies with complete 
mid- and far-IR spectroscopic coverage. This enabled 
us to transform the continuum luminosity functions into 
line luminosity functions up to z<4 and compute the line 
detection rates. 

2.1. Predicting continuum luminosity functions with the 
Franceschini et al. (2010) model 

As a first approach, we adopted the model developed 
by Franceschini et al. (2010). This is a backward 
evolution model that fits essentially all available data 
from Spitzer, ISO, COBE, and SCUBA. Moreover, it in- 
cludes constraints from preliminary results of the Her- 
schel Space Observatory surveys (presented in the Sci- 
ence Demonstration Phase papers of Oliver et al. 2010; 
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Fig. 1. — IR (8— 1000/im) continuum luminosity function con- 
structed from the Franceschini et al. (2010) model. 



Berta et al. 2010; Glenn et al. 2010; Gruppioni et al. 
2010; Nguyen et al. 2010). 

The model accounts separately for normal spiral galax- 
ies, actively star-forming galaxies, and for type 1 and 2 
AGNs. The actively star-forming population is further 
split into two galaxy classes, i.e., moderate- luminosity 
luminous infrared galaxies (LIRGs) and high-luminosity 
ULIRGs. The LIRGs have typical IR luminosities ~ 
10" Lq and the ULIRGs ~ lO^^ Lq. They are treated 
as two galaxy classes, with different luminosity functions, 
evolution rates, and spatial clustering properties. 

The broad-line AGNs-type 1 Seyferts and quasars- 
were modelled adopting luminosity functions and evolu- 
tion rates consistent with those observed in optical and 
mid-IR surveys (e.g. Rush, Malkan & Spinoglio 1993; 
Spinoglio et al. 1995). As detailed in Franceschini et al. 
(2010), important constraints on the infrared evolution 
properties of type 1 AGNs have been inferred from a flux- 
limited sample of 24/iTO-selected sources with complete 
spectroscopic classification, as reported by Rodighiero et 
al. (2010). These objects are easily identified by their 
flat spectral shapes over the optical through IR wave- 
length range (Spinoglio, Andreani & Malkan 2002). 

Narrow-line type 2 AGNs are instead much more dif- 
ficult to disentangle from starburst galaxies. Therefore 
their statistical properties and incidence among the IR 
population at high redshifts are still essentially unknown. 
At the present stage of knowledge, the model simply con- 
siders type 2 AGNs as a fraction of both low-luminosity 
and high-luminosity starbursts. In an attempt to con- 
strain such fractions, a follow-up analysis by Frances- 
chini et al. (2011, in preparation) has compared the star- 
formation and stellar-mass assembly histories of galaxies, 
and found that the two are consistent with each other if 
10% of the LIRG and 30% of the ULIRG objects are 
dominated by obscured (type 2) AGN accretion. These 



are the fractional contributions for type 2 AGNs adopted 
in the present work. It is clear from the above that this 
modelling might easily give results inconsistent with the 
value derived from both predictions by the AGN unifi- 
cation scheme and observations (e.g., from the SLOAN 
and FIRST surveys, Lu et al. 2010) of a --1 to 3 ratio of 
type-1 to type 2 objects. Indeed, the validity of the AGN 
unification at high redshifts has never been proven in de- 
tail, and will be one of the major outcomes of futures 
spectroscopic surveys in the IR. 

Figure 1 shows the continuum far-IR luminosity func- 
tions for all populations in each z range considered: 
0.<z<0.75, 0.75<z<1.25, 1.25<z<1.75, 1.75<z<2.25, 
2.25<z<2.75 and 2.75<z<4.0. 
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Fig. 2. — IR luminosity function from Gruppioni et al. (2011) 
for various galaxy populations, i.e., normal and starburst galaxies, 
low luminosity AGN, type 1 AGN and type 2 AGN. 



2.2. Predicting continuum luminosity functions with the 
Gruppioni et al. (2011) model 

Gruppioni et al. (2011) have also developed a back- 
ward evolution model fitting the main constraints pro- 
vided by IR/sub-mm surveys in the 15/im to 500/im 
range. In the mid-IR (MIR), it uses data from both 
ISO and Spitzer that are available in the literature, at 
I5fim from the ELAIS-Sl (Gruppioni et al. 2002), 
HDF-N, HDF-S and Marano fields (Elbaz et al. 1999), 
ultradeep lensed (Metcalfe et al. 2003), Lockman Deep 
and Shallow (Rodighiero et al. 2004), and at 24/im 
from the GOODS (Papovich et al. 2004) and SWIRE 
surveys (Shupe et al. 2008). In the FIR, we use data 
from Herschel, including those of the very recent FAGS 
Evolutionary Probe Survey (Berta et al. 2010; Grup- 
pioni et al. 2010), the Herschel-ATLAS Survey (Eales 
et al. 2010; Clements et al. 2010), and the Herschel 
Multi-tiered Extra-galactic Survey (Oliver et al. 2010; 
Vaccari et al. 2010). The model uses the classical ap- 
proach of evolving a local luminosity function in luminos- 
ity and/or density with z, with a different evolutionary 
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description for galaxies and AGN. What differentiates 
this model from others is the spectral energy distribu- 
tion (SED) scheme that it uses to distinguish between 
the different IR-bright galaxy populations. The scheme 
is based on a large spectroscopic study of MIR-selected 
sources (Gruppioni et al. 2008), and it divides sources 
into five (instead of four) broad SED classes. These are 
the normal spiral galaxies, the starburst galaxies, the ob- 
scured type 2 AGNs, the unobscured type 1 AGNs, and 
the objects containing a low luminosity (LL) AGN. In 
other words, the authors of this model have put partic- 
ular emphasis on determining the AGN contribution at 
low luminosities at all redshifts. The contribution of the 
various galaxy populations to the IR luminosity function 
can be seen in Figure 2. 




Fig. 3. — IR luminosity function from Valiante et al. (2009). 



2.3. Predicting continuum luminosity functions with the 
Valiante et al. (2009) model 

The third model that we used is the backward evolu- 
tion model of Vahante et al. (2009). It was developed 
using Spitzer and SCUBA observations, and it has been 
very successful in predicting Herschel results (Berta et 
al. 2010; Ahieri et al. 2010; Glenn et al. 2010; Oliver 
et al. 2010). This model allows us to take into account 
galaxies that are not 'pure' starburst or 'pure' AGN, and 
for which the ratio between IR lines might not be the 
ones expected assuming 'pure' SEDs. This is because 
the model considers all infrared galaxies as a single pop- 
ulation, assuming that starbursts and AGN coexist. It 
then uses an empirical relation to assign to each galaxy 
the fraction of the IR luminosity that is powered by the 
AGN for its given luminosity and redshift. This rela- 
tion was derived using a complete sample of local IRAS 
galaxies, and extrapolated to high z using Spitzer and 
SCUBA observations. Figure 3 shows the IR luminosity 
function predicted by this model for each z bin. 
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Fig. 4. — Critical density for coUisional de-excitation vs. ioniza- 
tion potential of IR fine-structure lines. 



2.4. Converting continuum luminosity Junctions to line 
or feature luminosity functions. 

The galaxy number counts per redshift and bolomet- 
ric IR luminosity bin that are predicted by each model 
need to be converted into line luminosity functions, in 
order to estimate the number of objects that will be 
detectable in various lines, and to assess whether sev- 
eral of the open questions presented in Section 1 can 
be addressed. For this purpose, we derived correla- 
tions between line and continuum luminosities, includ- 
ing lines for which such correlations were not previ- 
ously available in the literature, e.g., that of [Sill] for 
both AGN and star-forming galaxies. We examined the 
PAH feature at 11.25/im, the purely rotational H2 (0- 
0)S1 line at IT.OS/xm, and the [Nell] 12.8Aim, [NeV] 
14.3Aim, [Nelll] IS.S^im, [SIII] 18.7//m, [NeV] 24.3^m, 
[OIV] 26^1X1, [Sill] M.^,^lm, [OIII] 52^m, [NIII] ST^im, 
[01] ea^im, [OIII] 88/im, [Nil] 121.90Aim, [01] 145.52^m, 
and [CII] 157. 74^m fine-structure lines. These lines cover 
a wide parameter space of the critical density vs. ioniza- 
tion potential diagram (see Fig. 4), tracing different as- 
trophysical conditions: from photodissociation regions, 
to stellar/HII regions, to the AGN and coronal line re- 
gions (Spinoglio & Malkan 1992). This makes the com- 
bination of their ratios useful for the creation of AGN 
vs star-formation diagnostic diagrams (e.g. Spinoglio & 
Malkan 1992; Genzel et al. 1998; Dale et al. 2006; 
Smith et al. 2007). 

For lines at wavelengths shorter than 35/Ltm, we used 
the complete, 12/im-selected sample of local Seyfert 
galaxies (Tommasin et al. 2008, 2010) and the Bernard- 
Salas et al. (2009) sample of starburst galaxies to cali- 
brate the line luminosities to Lir . These samples have 
been extensively observed in the MIR with the IRS spec- 
trometer (Houck et al. 2004) onboard Spitzer (Werner et 



al. 2004), and the Spitzer spectra have been reduced and 
analysed in a consistent way. For the starburst galaxies, 
we excluded all objects for which there was evidence for 
the presence of an AGN from the literature or from the 
detection of [NeV] (see Table 1 of Bernard-Salas et al. 
2009). For the long- wavelength lines, we used the het- 
erogeneous sample of local galaxies compiled by Brauher 
et al. (2008) containing all observations collected by 
the LWS spectrometer (Clegg et al. 1996) onboard ISO 
(Kessler et al. 1996). The IR luminosities of the galaxies 
of our sample have been computed from the IRAS fluxes, 
using the formula of Lj j^^ representing the total mid- and 
far- infrared luminosity (Sanders & Mirabel 1996). All 
luminosities are in units of 10'''^ erg s~^. 

Using least-squares fitting, we obtained the following 
relations for the Seyfert galaxies, 
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relations are: 
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Iog(LpAH11.25) = 

with B. = 0.87, n 

log(i[NoII112.8l) = 

with R = 0.89, n 

Iog(L[NeV]14.32) = 
with R = 0.81, n 
log(L[NeIII]15.55) 

with R = 0.82, n 

log(L(H2)17.03) = 
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log(L[sillJ 18.71 ) = 

with R = 0.85, n 

log(^'[NoV]24.3l) = 

with R = 0.84, n 
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log(i[SIII]33.48) = 

with R = 0.91, n 

log{-^'[SiII]34.82) = 

with R = 0.89, n 



(0.95 ± 0.07) logiLin) - 

= 69, x"" = 7.1 
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= 87, = 8.8 
= (0.94 ± 0.08) Iog(L7H) - 

= 81, = 14.8 
= (0.95 ± 0.07) log(Ljfl) 

= 87, = 15-0 
(0.97 ± 0.05) log(LjH) - 
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(0.90 ±0.07) log(Liii) 
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(0.98 ± 0.08) log(Ljfl) 

= 71, x^ = 10.5 

= (0.88 ± 0.08) log{LiR) 
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(0.98 ± 0.05) log(LjH) 

= 75, x^ = 5.7 

(1.03 ± 0.06) log(L/H) 



(2.60 ±0.21) (1) 

-(3.25 ±0.18) (2) 

- (3.43 ± 0.25) (3) 

- (3.18 ±0.24) (4) 
-(3.79 ±0.15) (5) 
-(3.30 ±0.21) (6) 

- (3.45 ± 0.24) (7) 
-(2.66 ±0.25) (8) 
-(3.21 ±0.17) (9) 
-(3.14 ±0.20) (10) 
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where for each relation the Pearson coefficient R, the 
number of considered objects n and the computed 
are given. For the starburst galaxies, the corresponding 

® -L/jj is computed by fitting a siiiglc-tcmpcraturc dust cmis- 
sivity model (e oc u^^) to the flux in aU four IRAS bands, and 
should be accurate to ±5% for dust temperatures in the range 25 
- 65 K. We notice that the IR luminosities, as deflncd above, are 
model-dependent, and therefore could introduce some systemat- 
ics. However these do not affect the derived relations, as they are 
within the given errors. 



(1.17 ±0.11) log(Ljfl)- 
with R = 0.95, n = 14, x^ = 1-4 
log(L[Nen]i2.8i) = (1-17 ± 0.14) log(Ljfl) 
with R = 0.93, n = 14, x^ = 2.1 
log(L[Nem]16.B5) = (1.33 ± 0.18) log(LjH) 
with R = 0.90, n = 15, x^ = 3.9 
log(L(H2)i7.03) = (1-28 ± 0.14) logiLiu) - 
with R = 0.93, n = 15, x^ = 2.6 
Iog(i[Siiilis.7i) = (1-09 ± 0.15) log(LjH) - 
with R = 0.89, n = 15, x^ = 3.0 
log(L[oiv]25.89) = (1.24 ± 0.24) log(LjH) - 
with R = 0.85, n = 12, x^ = 1-6 
log(L[siII133.48) = (1-09 ± 0.10) log(LjH) - 
with R = 0.95, n = 15, x^ = 1-2 

log(i[SiII]34.82) = (1-11 ± 0.09) log(LjR) - 

with R = 0.96, n = 15, x^ = 0.92 

The best-fit solution is shown in Figure 
populations. Considering the sum of the 
derived the following generic relations: 

log(LpAHii.2.5) = (0.98 ± 0.06) log(LjH) - 
with R = 0.88, n = 83 x^ = 9.4 
log(L[NeII] 12.81 ) = (0.99 ± 0.06) log(LjK) - 
with R = 0.84, n = 101, x^ = 17.8 

log(L[NeIII]15.55) = (1-10 ± 0.07) log(Ljfl) 

with R = 0.83, n = 102, x^ = 24.2 
log(L(H2)i7.03) = (1.07 ± 0.05) log(Ljfl) - 



with R = 0.91, n = 91, x 
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log(i[siii]i8.7i) = (0.97 ± 0.06) log(Ljfl) - 
with R = 0.88, n = 70, x^ = 7.9 
log(i[OIV125.89) = (0.95 ± 0.11) log(LjH) ■ 
with R = 0.68, n = 95, x^ = 36.6 
log(L[siiii33.48) = (0.99 ± 0.05) log(LrH) - 
with R = 0.92, n = 90, x^ = 7.4 
log(L[SiII134.82) = (1.04 ± 0.05) log(LjH) - 
with R = 0.91, n = 87, x^ = 8.4 



(3.03 ±0.32) (11) 

- (3.65 ± 0.40) (12) 

- (4.85 ± 0.52)(13) 
(5.10 ±0.42) (14) 

- (3.79 ± 0.45) (15) 

- (5.13 ±0.74) (16) 

- (3.35 ± 0.29) (17) 

■ (3.26 ± 0.25) (18) 

5 for each of the 
populations, we 

(2.69 ±0.18) (19) 

- (3.26 ± 0.20) (20) 

- (3.72±0.23)(21) 
(4.19 ±0.16) (22) 
(3.47 ±0.20) (23) 

■ (3.04 ± 0.34) (24) 
(3.21 ±0.14) (25) 
(3.15 ±0.16) (26) 



For the far-IR lines (Fig. 6), we obtain: 

log(i[Oiiil5i.8i) = (0.88 ± 0.10) log(L/fl) - (2.54 ± 0.31) (27) 
with R = 0.91, n = 16, x^ = 2.6 

log(L[Niii]57.32) = (0.78 ± 0.10) log{Lin) - (2.58 ± 0.32) (28) 
with R = 0.94, n = 10, x^ = 0.14 

log(L[oi]63.i8) = (0.98 ± 0.03) log(Ljfl) - (2.70 ± 0.10) (29) 
with R = 0.94, n = 109, x^ = 9-1 

Iog(L[oiii]88.36) = (0.98 ± 0.10) log(Ljfl) - (2.86 ± 0.30) (30) 

with R = 0.81, n = 55, x^ = 12.7 

log(L[Nii] 121.9) = (1.01 ± 0.04) log(L/fl) - (3.54 ± 0.11) (31) 
with R = 0.93, n = 100, x^ = 13.3 

Iog(i:[oi]i45.5) = (0.89 ± 0.06) log(Ljfl) - (3.55 ± 0.17) (32) 
with R = 0.91, n = 46, x^ = 10-0 

log(L[cii] 167.7) = (0.89 ± 0.03) log(LjH) - (2.44 ± 0.07) (33) 
with R = 0.92, n = 217. x^ = 42.3 
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log[L(FIR)/10<'erg s-'] log [L(FIR)/10<'erg s-'] 



Fig. 5. — Correlations between the various feature and line luminosities and the far-IR luminosity for the Seyfert galaxies of the complete 
12/im galaxy sample (Tommasin et al. 2008, 2010), and for the pure starburst galaxies of the sample of Bernard-Salas et al. (2009). 
The dotted, broken and solid lines represent the least-squares fit of the data of the Seyfert, the pure starburst galaxies and all galaxies 
populations together, respectively. Figures (c) and (g) have only the Seyfert galaxies, because the [NeV] lines are not detected in starburst 
galaxies. 



FIR spectroscopic cosmological surveys 7 
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Fig. 6.— Correlations between the [OIII]52/im, [NIII]57/im, [OI]63/im, [OIII]88/im, [NII]122/im, [OI]145Atm and [CII]158/im luminosity 
and the far-IR luminosity for the galaxies observed with the ISO-LWS spectrometer (Brauher et al. 2008). 



For all correlations the hypothesis that the variables 
are unrelated can be rejected at a level of significance 
which is always less that 10~^. Wu et al. (2010) pre- 
sented the relation between the total infrared luminosity 
and the PAH emission band at 11.25/im for AGN and 
starburst galaxies of the 24/im flux limited intermediate 
redshift (< z >~ 0.14) sample of 5MUSES (Helou et 
al. 2011, in preparation). Their result is comparable 
to ours. By inverting our relations we derive a slope of 
1.05±011 for the Seyfert galaxies and 0.85±0.08 for the 
starburst galaxies, compared to their slopes of 1.00±0.04 
and 0.98±0.03, respectively. 



3. RESULTS 

3.1. Number counts per spectral line 

To compare and visualize our results for the three 
galaxy-evolution models, we need to adopt a line- 
detection sensitivity curve as a function of wavelength A, 
an integration time, and a field of view for the simulated 
observations. For this purpose, we opt to use numbers 
relevant to future missions or facilities. For the FIR do- 
main, we use the sensitivity curve proposed for SPICA's 
far infrared instrument (SAFARI), while for the submil- 
limeter domain, we use the sensitivity curve that has 
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been estimated for a R=1000 resolution spectrometer at 
the focal plane of the CCAT telescope. Details on these 
instruments are presented in the Appendix. Moreover, 
we selected a common integration time of 1 hour and a 
total field of view to be covered by our simulated survey 
of 0.5 deg^. For an instrument, such as SAFARI, with a 
2' X 2' field of view, this corresponds to 450 hours of in- 
tegration time, to be compared to 4.5 hours for a CCAT 
spectrometer, assuming a field of view of 20' x 20'. 

The resulting number of AGN and starbursts that will 
be detectable in each line as a function of z with SA- 
FARI is presented in Tables 1 and 2, respectively, for 
the Franceschini et al. (2010) model. The same results 
are presented in Tables 3 and 4 for the Gruppioni et al. 
(2011) model. The number of detectable galaxies based 
on the Valiante et al. (2009) model is presented in Ta- 
ble 5. We note that for the sake of completeness, and 
to assist further planning and designing of new instru- 
mentation, we also present in Tables 1-5 the predicted 
number of sources that are detectable in the low— Z bins 
for the short-wavelength lines and those in the high-Z 
bins for the long-wavelength lines, even lines outside the 
nominal SAFARI spectral range. We have simply as- 
sumed a fiat extrapolation of the SAFARI sensitivities 
to shorter and longer wavelengths. Tables 6, 7, and 8 
present the simulation results for CCAT. 

A basic result of this analysis is that the total number 
of detectable objects agrees to within a factor of 2—3 for 
most lines and z ranges, and that at least a thousand 
galaxies will be simultaneously detected in four lines at 
5(7 over a half square degree. A comparison of the output 
of the three models is plotted in Fig. 7 for SAFARI. A 
survey of the assumed sensitivity will detect bright lines 
(e.g., [O I] and [O III]) and PAH features in thousands 
of galaxies at z>\. Hundreds of 2;>1 AGN will be de- 
tected in the [O IV] line, and several tens of z>\ sources 
will be detected in [Ne V] and H2. For H2 in particular, 
this number corresponds to a lower limit. Our models do 
not account for an increase of the H2 emission efficiency 
as a cooling mechanism, or for the H2 mass content 
to increase with increasing z and decreasing metallic- 
ity. Neither do the applied continuum-to-line luminosity 
relations include sources of extremely high L(H2)/Lir 
ratios, associated with shock fronts due to galaxy colli- 
sions or AGN feedback mechanisms (Cluver et al. 2010; 
Ogle et al. 2010). 



Our line detectability results are sensible when com- 
pared with predictions from local galaxy templates. 
To make this comparison, we used four objects with 
well-determined MIR and FIR spectra. These are 
NGC1068, a prototypical Seyfert 2 AGN, NGC4151, a 
well-studied type 1 AGN, the prototypical moderate- 
luminosity starburst M82, and a starburst-dominated 
ULIRG, IRAS17208-0014. Their line intensities are 
taken from Alexander et al. (2000); Spinoglio et al. 
(2005), Sturm et al. (1999); Spinoglio et al. (1997), 
Farrahetal. (2007); Brauher et al. (2008) and Forster 
Schreiber et al. (2001); Colbert et al. (1999), respec- 
tively for the four templates. We scaled the bolometric 
IR luminosity of all systems to 10^^ L0, and show which 
lines can be observed as a function of z in Figs. 8 and 
9. This basic comparison confirms that faint lines can be 
detected in z^\ ULIRGs with SAFARI. Emission from 
[O IV] wfil be observable in z-3 ULIRGs with SAFARI, 
while CCAT wifi observe [C II] out to z-5. 

In Fig. 10, and 11 we further investigate in what type of 
IR-bright galaxies can each line be detected as a function 
of z. We only make this comparison for the Franceschini 
et al. (2010) model, as the three models do not differ 
in terms of their IR-luminosity classification, and as the 
total number of sources is comparable in all cases. We 
find that for the brightest lines, such as [Si II], SAFARI 
will be able to observe LIRGS even at z>3. However, for 
the typical [Ne II] and [O III] 52 /xm line luminosities, 
the transition from LIRGs to ULIRGs will occur at z~2, 
and for the [Ne V] and H2 SI lines at z^\. CCAT wiU 
be highly complementary to SPICA, as it will be able to 
observe the [O III] 88 /im line at z>1.3, where this line 
leaves the spectral range of SAFARI. We also find that 
CCAT will be a most efficient instrument for studies of 
[C II] , an important coolant of the interstellar medium, 
at all z<5. At 3<z<4 alone, it will detect more than 300 
galaxies at ha level in a 0.5 deg'^ survey. 

3.2. Line luminosity function predictions in the 
IR/suhmm 

We present in Fig. 12 and 13 the predicted luminos- 
ity functions of AGN and starburst galaxies for each line 
and feature in the Franceschini et al. (2010) model. It is 
clear from the figure that the AGN line luminosity func- 
tions (dashed lines in the figure) for the lowest rcdshift 
range (0.<z<0.75) are in agreement with the local lumi- 
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Fig. 8. — a) Predictions of lines observable with SPICA/SAFARI, based on local templates, scaled to an intrinsic luminosity of 10^^ Lq. 
Selected diagnostic lines are shown as a function of z for a type 1 AGN (NGC4151) and a type 2 AGN (NGC1068). For comparison, we 
overplot the 5cr (in 1 hour) sensitivity threshold of SPICA-SAFARI as a function of wavelength (see Appendix). For completeness, we also 
report the sensitivities for Herschel PACS (Poglitsch et al. 2010), for the Mid-Infrared Instrument (MIRI, Wright et al. 2004), which will 
be onboard the James Webb Space Telescope (Gardner et al. 2006) and which will operate from 5 to 28^m, and for the Mid-ir Camera and 
Spectrograph (MCS) planned for SPICA (Wada & Hirokazu 2010). b) Same as panel a) for a moderate-luminosity prototypical starburst 
(M82), and a starburst-dominated ULIRG (IRAS17208-0014). 




Fig. 9. — Same as Fig. 8, for long-wavelength lines, where the complementarity between SPICA and CCAT is clearly shown. For the 
adopted sensitivities of CCAT we refer to the Appendix. 
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Fig. 10. — Prediction of the number of sources of a 0.5 deg^ spectroscopic survey with SAFARI based on Franceschini et al. (2010), 
giving the number of detectable starburst galaxies (divided by Lir) and AGN (divided by obscuration) at the 3(t level. The adopted line 
flux sensitivities as a function of wavelength are given in the Appendix. The left panels correspond to the AGN predictions, the right 
panels to the starburst predictions. As for the former, the number of type 2 AGNs associated with the LIRG and the ULIRG populations 
are shown separately (in black and red, respectively). 



nosity functions, at an average < z > of 0.03. The space 
densities of starburst galaxies are expected to be higher 
than those of AGN, for any hue, except for the [OIV] 
hue, which~as is weh known-is much fainter in starburst 
galaxies. The volume densities of AGN drop faster with 
z than those of starbursts when traced by the [Ne II] and 
[Ne III] lines, possibly due to a saturation of the lines in 
AGN of high luminosities. The total number of AGN 
detected in the same lines is one to two orders of mag- 



nitude lower than that detected for starbursts at any z. 
The number of sources detected in [O III] 88 /xm is com- 
parable for AGN and starbursts at all z. This result most 
likely differs from that for [Ne III] which comes from an 
ion of comparable ionization potential to O III, because 
of the single relationship used to convert the line luminos- 
ity to iiR, for long wavelength lines. This is even true 
for the [O I] 63 /im line, which also has a high critical 
density for collisional de-excitation, ~10^ cm~^, mak- 
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Fig. 11. — Prediction of the number of sources of a 0.5 deg^ 
spectroscopic survey based on Franceschini et al. (2010), giving 
the number of detectable galaxies at the Scr level with CCAT. The 
adopted line flux sensitivities as a function of wavelength are re- 
ported in the Appendix. 
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Fig. 12. — Predicted line luminosity functions of [NeV]14.3/.tm, 
[NeIII]15.5AJm, [NcV]24.3AJm, [OIV]26AJm, [OIII]52Atm and 
[OIII]88/im, for SAFARI. Dashed lines correspond to the predic- 
tions for AGN, while solid lines correspond to the predictions for 
starbursts. Where available, the comparison with the observed 
local LF of AGNs from Tommasin et al. (2010) is given. 



ing it bright in AGN. The contamination of the [O IV] 
AGN luminosity functions from starbursts is minimal all 
the way through z—4:. The volume density of actively 
accreting black holes as traced by [O IV] and [Ne V] in- 




L(line) {erg s' 



Fig. 13. — Predicted line luminosity functions of PAH 11.25^im, 
[NeII]12.8Atm, H2 17/im, [SiII]34.8Atm , [NIII]57Atm and [Ol]63^lra, 
for SAFARI. Dashed lines correspond to the predictions for AGN, 
while solid lines correspond to the predictions for starbursts. 
Where available, the comparison with the observed local LF of 
AGNs from Tommasin et al. (2010) is given. 



creases up to i^z<3 before it drops back down at z—A. 
The same applies for tracers of star-formation, reproduc- 
ing the suggested coevolution of black hole growth and 
stellar mass build-up. 

4. DISCUSSION: NEW PARAMETER-SPACE COVERAGE 
BY SPICA AND CCAT 

Which of the questions raised in the Introduction 
will the proposed future telescopes address? The lines 
predominantly emitted by ions in star-forming com- 
plexes, like [Si II], [C II], and [Ne II], will be detected in 
Li^>10^^ Lq systems at least out to z~2. They will also 
be detected in 10^° Lq< Lir<10^^ Lq galaxies at least 
out to This result indicates that star- formation 

tracers will be detected and compared for the sources 
that are mainly responsible for the formation of present- 
day ellipticals. The creation of line-ratio diagnostic di- 
agrams, and the comparison of the line to IR contin- 
uum luminosities in hundreds of sources will help us, in 
combination with imaging data, to further address the 
star-formation bimodality, and to obtain a more coher- 
ent picture of the intermediate/high-z IR-bright popula- 
tions. Several sources will be detected with CCAT even 
past the peak of star-formation activity at z>3 in the 
strong [C II] line, which will help us constrain the shape 
of Lilly-Madau diagrams at such high rcdshifts. 

To date, the detection of FIR fine-structure lines has 
only been achieved in lensed z>l systems. Herschel 
PACS detected the [O III] 52^m line with a flux of 9x 
10~i^Wm"2 jj^Ag F10214+4724 at 2=2.28 (Sturm et 
al. 2010). In a z=1.32 source, MIPS JI42824.0+352619, 
[O III] 52/im and [O I] 63/im were detected with fluxes 
of 3.7 and 7.8x lO^^^Wm"^, respectively (Sturm et al. 
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2010) . The [O III] 88,um line was detected at the Cahech 
SubmiUinieter Observatory (CSO) with the ZEUS spee- 
trometer (Ferkinhoff et ah 2010b) in APM 08279+5255 
at z=3.9 and SMM J02399-0136 at ^=2.8 (Ferkinhoff et 
al. 2010a), with fluxes of 2.68 and 6.04 xlQ-i^Wm-^. 
All four systems are lensed galaxies, with magnification 
factors estimated to be in the range 2.4—90 (Egami et 
al. 2000; Ao et al. 2008; Riechers et al. 2009; Ivison 
et al. 2010). Such experiments were not possible for 
unlensed galaxies with the present-day missions, leaving 
ample room for new discoveries for SPICA and CCAT. 

The recent detection of [NII]122/im (Ferkinhoff et al. 

2011) and [CII]158Aim (Stacey et al. 2010) lines in a 
few high redshift galaxies with the ZEUS spectrometer 
(Ferkinhoff et al. 2010b) at the CSO shows that these 
lines can be much brighter than in local galaxies. The 
[NII]122/im line to FIR luminosity ratio is 2-10 times 
higher in the two observed galaxies (H1413+117 and 
SMMJ02399-0136) than in the local galaxies that we used 
for deriving the line to continuum relations adopted for 
our predictions. Similarly the [CIIJ/FIR luminosity ra- 
tio in the 1 < z < 2 galaxies can be twice as high as 
the value observed in the local Universe, despite hav- 
ing a large rage from 0.024% to 0.65% . This indicates 
that our estimates can be considered conservative, and 
that CCAT could observe many more galaxies than what 
shown in our predictions. It could potentially trace the 
star formation history of galaxies back to very high red- 
shift. 

A simple conversion of the SAFARI sensitivity limit 
into an H2 line flux and mass indicates that detection 
of H2-bright galaxies at z>6 will be feasible over large 
areas. An SI flux of 2.5x 10^^^ W m^^ at a typical tem- 
perature of 300K, yields an H2 mass of 3xlO^"Af0 at 
z=6. Cosmological simulations indicate that 2~6 galax- 
ies of such a molecular gas mass can exist (Obreschkow & 
Rawlings 2009). Unless there is a significant increase in 
the efficiency of H2 to cool the ISM, or unless there is a 
considerable number of sources with an AGN jet-ISM in- 
teraction that leads to very high L( H2)/ Xir ratios (Ogle 
et al. 2010), our conclusion will not hold at ^;~10. At 
that z, the same temperature and SI flux correspond to 
a mass of 10" Af©. Even though the H2 (O-O)SO line 
at 28.03 /im will be in the spectral range of CCAT at 
z> 6, its detection will be equally hard: assuming an SO 
flux of 2x 10~^^ W m~^ and an ortho-to-para ratio of 3, 
the minimum detectable H2 mass would be ^10"'^'^ Mq at 
z=6. SPICA will thus be unique for H2 studies seeking 
the end of the reionization era. 

Another parameter space unique to SPICA will be 
the detection of resolved [O IV] 25.89 jum emission up 
to 2;~4. The calibration of the widths of high-ionization 
MIR lines, as probes of the narrow-line-region kinemat- 
ics, to the black hole mass was attempted with Spitzer 
for sources out to z=0.3 (Dasyra et al. 2008, 2011). 
It might be useful for the mid-infrared instrument of 
JWST, which will be able to observe the [Ne V] 14.32 /xm 
line out to 2;~0.8. The breakthrough for such studies will 
come however from SPICA. SAFARI will access the z>l 
Universe and go past its accretion-rate-history peak. The 
shape of the accretion rate functions will be then com- 
pared to the shape of the mass functions of the obscured 
and unobscured black holes upon which the material is 
accreted. The inclusion of obscured black holes in mass 



functions of unobscured AGN (e.g., Vestergaard & Os- 
mer 2009) could change their observed volume-density 
normalization factors and shapes. The most obscured 
AGN are thought to be missing from a redshift of only 1 
(e.g., Gilh et al. 2007), unlike the bulk of the black hole 
growth that is thought to be occurred at 2; > 1. 

5. SUMMARY 

We used three galaxy formation and evolution mod- 
els, constrained by luminosity functions from multi- 
wavelength observations, to predict the number of star- 
forming and active galaxies that are observable in any IR 
luminosity bin with look-back time. We converted the 
number counts of galaxies per IR luminosity to number 
counts of galaxies per line luminosity using several (new) 
line-to-continuum conversion relations, built upon local 
AGN and starburst galaxy samples. We compute our re- 
sults for an hour-long integration/FoV half-square-degree 
survey, and for the sensitivity values of SPICA SAFARI 
and CCAT. Their anticipated values are 2.5x10"^^ W 
m-2 at 160 ^m and 1.1x10"^^ W m'^ at 620 /xm, re- 
spectively. These telescopes/instruments were selected 
because they are designed to be survey machines, able 
to perform large blind cosmological surveys in reason- 
able integration times. We find that SAFARI will de- 
tect thousands of 2;>1 galaxies in bright low-ionization 
fine-structure lines such as [Sill] and [Nell], and sev- 
eral tens of 1 < 2; < 2 galaxies in fainter lines such as 
[O IV], [Ne V] and H2 (O-O)Sl. For the bright fines, nor- 
mal galaxies will be observed out to z^l, LIRGs out 
to 2;~2, and ULIRGs to even higher z. This means 
that studies of the ionized gas properties in the galax- 
ies that form the present day massive ellipticals will be 
feasible. AGN/star- formation diagnostic diagrams will 
be obtained for different classes of IR-bright galaxies at 
1 < z < 2, which will enable us to not only look for 
a redshift evolution of line ratios, but also for a lumi- 
nosity evolution within each z range. Further tests of 
the black hole growth - galaxy build-up coevolution sce- 
nario will be performed, as the creation of accretion-rate 
functions and mass functions will be determined out to 
2:^4, for both obscured and unobscured black holes using 
the [O IV] line. Over large areas, SAFARI could also be 
able to detect H2 at z>6, and help constrain the end of 
the reionization era. In the light of the current findings 
at ground-based submillimeter telescopes of substantially 
brighter fine structure lines in high redshift galaxies com- 
pared to local galaxies, we are even more confident that 
CCAT will be unique for studying the star formation his- 
tory of galaxies back to very high redshift. 
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APPENDIX 

DESCRIPTION OF THE SPICA MISSION AND THE CCAT FACILITY TO BE USED FOR FUTURE, 
BLIND COSMOLOGICAL SURVEYS IN THE IR AND SUBMM 

The deep cosmological surveys undertaken by ISO (Kcssler ct al. 1996), Spitzer (Werner et al. 2004), AKARI 
(Murakami et al. 2007; Goto et al. 2010), WISE (Wright ct al. 2010), and 77ensfi/ie/ (Pilbratt ct al. 2010) wiU have 
produced catalogues containing the fluxes of many tens of thousands of IR-bright sources by the early 2020's. These 
catalogues will provide excellent targets to be followed up by ALMA (Brown, Wild & Cunningham 2004; Wootten 
2008), JWST (Gardner et al. 2006), SPICA (Swinyard et al. 2009), and CCAT (Sebring 2010). Among the listed 
facilities, JWST and ALMA will be more suited to deep follow-up spectroscopy of known targets. SPICA and CCAT 
will be suited for performing blind large-scale spectro-photometric surveys, because of the wide field of view (of several 
arcminutes squared) of their instruments. Covering different wavelength ranges, these two instruments will be highly 
complementary. In the rest of the Appendix, we provide details of their present-day design concepts. 

SPICA is a proposed JAXA-led astronomical mission with suggested contributions by European, Korean, and 
possibly US institutions, and with a launch date planned in the early 2020s (Swinyard et al. 2009). With a 3-m 
mirror that is actively cooled to <6K, and a state-of-the art focal plane instrument suite, SPICA will make imaging and 
spectroscopic observations over the 5-210/Ltm spectral range with unprecedented sensitivity. It will offer an improvement 
in raw photometric sensitivity with respect to Herschel of two orders of magnitude in the FIR. In the MIR, it will 
extend the capabilities of JWST with an uninterrupted spectral coverage in the range b-SSfim. The Mid Infrared 
Camera and Spectrometer (MCS) planned for SPICA (Wada k Hirokazu 2010; SPICA Study Team CoUaboration 
2010) is an integral field unit with a field of view of 12" x 6" at 10-20^m and 12" x 12.5" at 19.5-36.1/im. Its 
5 cr, 1 hour sensitivity is in the range 2-2.5 xlO^^'^Wm^^ (Wada & Hirokazu 2010). Observing capability in the 
FIR is provided for SPICA by SAFARI (Swinyard et al. 2009; SPICA Study Team Collaboration 2010). Proposed 
by a consortium of European institutes (with Canadian and Japanese participation), SAFARI is an imaging Fourier 
Transform Spectrometer (FTS) with a field of view of 2' x 2'. The FTS provides an instantaneous spectral coverage of 
the 34 - 210/im wavelength range and spectral resolution modes with A/AA of 2000 (at 100 fJ,m), ^ few hundred, or even 
as low as 20 < A/AA < 50. With its sensitive superconducting transition edge sensor (TES) detectors (Khosropanah et 
al. 2010), SAFARI will offer a factor of ~100 increase in raw sensitivity in the continuum, and ~15 in high-resolution 
mode spectroscopy compared to Herschel PACS (Poglitsch et al. 2010). The improvement in spectral mapping speed 
over PACS wiU be of more than 100 at A/AA~ 2000. Assuming a Noise Equivalent Power (NEP) of 2 x IQ-^^W/VHz 
for the SAFARI detectors, the 5 cr, 1 hour detection limits are predicted to be 4.16, 2.58, 1.89 and 2.48 xlO"^^Wm"^ 
for the four planned spectral bands centred at 48, 85, 135 and 160/Ltm, respectively, for an unresolved line at A/AA~ 
2000 (Swinyard, B., 2011. priv. comm.). 

The Cerro Chajnantor Atac;ama Telescope (CCAT) will be a 25m submm wave telescope to be constructed near 
the summit of Cerro Chajnantor (Sebring 2010). The CCAT science case, observatory requirements, and conceptual 
design were developed as part of a study jointly funded by Cornell and Caltech/JPL, which resulted in a Feasibil- 
ity/Concept Design Study Report, available at: http://www.submm.org/doc/2006-01-ccat-feasibility.pdf. The design 
of the CCAT telescope is currently undergoing changes to increase the field of view from 20 x 20 arcminutes square to 
1 square degree (Sebring 2010). This would substantially increase the survey capability of the telescope. CCAT will 
carry out spectroscopic surveys of submm galaxies, using multi-object versions of broadband direct-detection grating 
spectrometers such as Z-Spec (Bradford et al. 2004) and ZEUS (Ferkinhoff et al. 2010b), now in use at the CSO. 
Conceptual development indicates that spectrometers capable of observing 10-100 objects simultaneously while span- 
ning multiple atmospheric windows will be feasible (Stacey et al. 2006) . Sensitivity estimates based on the 25 m 
telescope CCAT telescope with a 10 /im rms surface on Cerro Chajnantor (5600 m elevation) for a spectrometer with 
a resolution of 1000 (5(j, 1 hour) that include wavelength-dependent, typical precipitable water vapor corrections are: 
2.0, 1.9 and 1.0 xlO-^^Wm-^ at 200Mm, 230/im and 291/im respectively, 1.8, 1.3 and 1.1 xlO^i^Wm"^ at 350Atm, 
450/im and 620/im respectively and 3.8 and 2.6 xl0~^°Wm~^ at 740/im and 865/im respectively (Stacey, G. 2011, 
private comm.). 
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TABLE 1 

Number of AGN detectable in a SAFARI survey of 0.5^ in IR lines as a function of redshift at 5cr (3ct) in 1 hr. 

INTEGRATION PER Fo\'. EOLLOWIXC FRAXCE^rniNI ET AL. ^20101 



line/redshift 


0<z<0.75 


0.75<z<1.25 


1.25<z<1.75 


1.75<z<2.25 


2.25<z<2.75 


2.75<z<4 


all redshifts # 


PAH(11.25jum) 




^;qq ('\ a An \ 4- 
ooo. ^^i44y._) y 


A'in /^'7nQ'\ + 

4oy. \iyio) \ 


04.U ^^iott.j 


K QC ( A'7 '7\ 

O.ttO \^:( .i ) 


■ ■ ■ (1.35) 


OU.U (^icS/.J if 




Nell] 12.81//m 


152. (366.) t 


122. (247.) t 


14.4 (45.0) 


90 (6 75) 




... ( ... ■) 


15.3 (52.2) X 




NeV] 14.32 /im 


43.2 (152.) t 


20.2 (82.3) t 


0.45 (5.85) 


... ( ... ) 


... ( ... ) 


... ( ... ) 


0.45 (5.85) t 




Nelll] 15.55/im 


152. (366.) t 


35.1 (82.3) 


5.85 (45.0) 


0.45 (1.80) 


■ • ■ (0.45) 


... ( ... ) 


41.4 (130.) t 


H2(17.03/im) 


10.8 (43.2) t 


■■■ (6.75) 


■ ■ ■ (0.45) 


... ( ... ) 


... ( ... ) 


... ( ... ) 


• ■ ■ (7.20) t 




siii] ls.^l^lm 


43.2 (106.) t 


3.15 (11.2) 


■ ■ ■ (0.90) 


... ( ... ) 


... ( ... ) 


... ( ... ) 


46.3 (118. )t 




NeV] 24.32^4111 


26.6 (69.8) 


11.2 (55.8) 


14.4 (45.0) 


0.90 (6.75) 


■ ■ ■ (0.90) 


... ( ... ) 


53.1 (178.) 




OIV] 25.89Atm 


152. (366.) 


82.3 (176.) 


73.8 (246.) 


17.1 (54.0) 


0.90 (14.8) 


... ( ... ) 


326. (857.) 




SIII] 33.48Atm 


69.8 (210.) 


122. (247.) 


45.0 (174.) 


6.75 (32.8) 


0.90 (14.8) 


... ( ... ) 


244. (679.) 




Sill] 34.8lAtm 


210. (366.) 


333. (633.) 


174. (439.) 


54.0 (204.) 


14.8 (121.) 


1.35 (11.7) 


787. (1775.) 




OIII] 51.81^111 


464. (687.) 


333. (633.) 


246. (563.) 


54.0 (204.) 


14.8 (28.8) 


• ■ ■ (0.90) t 


892. (2117.) t 




NIII] 57.32^tm 


106. (281.) 


55.8 (176.) 


14.4 (73.8) 


• ■ ■ (0.90) 


... ( ... ) 


... ( ... ) 


176. (532.) 




OI] 63.18/im 


687. (945.) 


1128. (1817.) 


563. (863.) 


204. (381.) 


47.7 (179.) 


4.50 (22.5) t 


2833. (4184.) J 




OIII] 88.35;im 


687. (945.) 


459. (856.) 


246. (563.) 


88.2 (204.) t 


28.8 (77.4) t 


0.90 (11.7) t 


1392. (2364.) j 



Notes: total number of AGN inside detectable the SAFARI spectral range; t- 
outside the SAFARI spectral range. 



outside the SAFARI spectral range;J: excluding detections 



TABLE 2 

Number of starburst galaxies detectable in a SAFARI slirvey of 0.5^ in IR lines as ,\ function of redshift at 5<t (3a) in 

1 HR. integration PER FoV, FOLLOWING FrANCESCHINI ET AL. (2010) 



linc/redshift 



0<z<0.75 



0.75<z<1.25 1.25<z<1.75 1.75<z<2.25 2.25<z<2.75 



2.75<z<4 



all redhifts # 



PAH(11.25/im) 
[Nell] 12.81/im 
[Nelll] 15.55/im 
H2(17.03^/m) 
S'ln] 18.7lAtm 
OIV] 25.89/im 
5777] 33.48/^m 
5i77] 34.81/im 
0777] 51.8lMm 
Af777] 57.32/im 
OI] 63.18/im 
0777] 88.35//m 



1023. (1519.) t 
636. (1023.) t 
115. (254.) t 
28.8 (43.2) t 
176. (355.) t 
2.70 (9.45) 
355. (636.) 
(816.) 
(1257.) 
(482.) 
(1808.) 
(1808.) 



482. 
816. 

176. 
1257. 
1257. 



974. (1790.) 1 
692. (1340.) 1 

80.5 (118.) 

9.90 (29.7) 

50.8 (118.) 

1.80 (9.90) 

495. (974.) 
974. (1790.) 

495. (974.) 

80.5 (258.) 
1790. (2960.) 
692. (1340.) 



608. (984.) t 
239. (460.) 
37.3 (101.) 
2.25 (7.65) 
19.3 (61.6) 
1.35 (7.65) 
337. (608.) 
608. (984.) 
337. (783.) 
19.3 (101.) 

783. (1221.) 
337. (783.) 



119. (277.) 
72.9 (187.) 
9.00 (44.1) 
0.45 (1.80) 
2.70 (22.9) 
0.45 (0.90) 
119. (277.) 
277. (519.) 
72.9 (277.) 
■■■ (1.80) 
277. (519.) 
119. (388.) t 



39.1 (105.) 

20.2 (64.3) 
1.35 (8.10) 
... ( ... ) 

1.35 (8.10) 
... ( ... ) 

64.3 (164.) 
164. (339.) 
2.25 (39.1) 
... ( ... ) 

64.3 (242.) 
39.1 (105.) t 



1.80 (15.7) 
0.90 (6.30) 
0.45 (1.80) 
... ( ... ) 

•■• (0.45) 
... ( ... ) 

6.30 (30.1) 

30.1 (81.4) 
• • ■ (0.90) t 
... ( ... ) 

6.30 (30.1) t 
0.90 (15.7) t 



160. (398.) t 
333. (718.) t 
129. (273.) t 
12.6 (39.1)t 
250. (566.)t 
6.30 (27.9) 
1377. (2689.) 
2535. (4529.) 
1723. (3330.) t 

276. (843.) 
4171. (6750.) t 
2286. (3931.) J 



Notes: #: total number of AGN inside detectable the SAFARI spectral range; f: outside the SAFARI spectral range;^: excluding detections 
outside the SAFARI spectral range. 

TABLE 3 

Number of AGN detectable in a SAFARI survey of 0.5^ in IR lines as a function of redshift at 5<t (3it) in 1 hr. 

INTEGRATION PER FoV, FOLLOWING GRUPPIONI ET AL. (2011) 



line/redshift 



0<z<0.75 



0.75<z<1.25 



1.25<z<1.75 1.75<z<2.25 2.25<z<2.75 



2.75<z<4 



all redhifts # 



PAH (11. 25^01) 
[Afe77] 12.8lAtm 
[NeV] 14.32 ^lUl 
[Afe777] 15.55^111 
H2 (17.03^111) 
'SIII] IS.nfim 
NeV] 24.32ftm 
OIV] 25.89Mm 
5777] 33.48Atm 
Sj77] 34.8lMm 
0777] 51.8lMm 
Af777] 57.32/im 
OI] 63.18^111 
0777] 88.35Aim 



1057.(1451.) t 
398. ( 701.) t 
177. ( 398.) t 
398. ( 701.) t 
82.3( 177.) t 
177.( 398.) t 

124. ( 240.) 

491. ( 701.) 

240.( 491.) 

491. ( 701.) 

814.(1184.) 

398. ( 701.) 
1057.(1316.) 
1057.(1316.) 



1076. (2317.) t 
111. ( 326.) t 
9.90 ( 58.0) t 
16.6 ( 58.0) 
0.45 ( 3.15) 
1.80 ( 5.40) 
5.40 ( 32.8) 
58.0 ( 192.) 
111. ( 326.) 
512. (1076.) 
512. (1076.) 
32.8 ( 192.) 
1862. (2805.) 
762. (1443.) 



486. (1126. )t 
9.90 (28.8 ) 
0.90 (5.85 ) 
5.85 (28.8 ) 
■ ■ ■ (0.90 ) 
0.45 (1.80 ) 
9.90 (28.8 ) 

47.7 (187. ) 

28.8 (121. ) 
121. (486. ) 
188. (759. ) 
9.90 (47.7 ) 
759. (1593.) 
188. (759. ) 



64.3 (148.) 
3.60 (12.6) 
• ■ ■ (0.90) 
1.80 (6.30) 
... ( ... ) 

■ • • (0.45) 
3.60 (12.6) 
22.0 (64.3) 
12.6 (39.1) 
64.3 (209.) 
64.3 (209.) 
0.45 (3.60) 
209. (391.) 
100. (209.)t 



13.9 (75.6) 
0.45 (1.80) 
... ( ... ) 

0.45 (1.80) 
... ( ... ) 

■■■ (0.45) 
0.45 (3.60) 
3.60 (25.2) 
6.75 (25.2) 
45.9 (170.) 
6.75 (45.9) 

■•• (0.45) 
75.6 (230.) 
45.9 (117.) t 



0.90 (7.65) 
( •■■ ) 
( •■■ ) 
(0.45) 

( ■■■ ) 
( ■■■ ) 
( •■■ ) 
0.90 (4.05) 
0.90 (4.05) 
7.65 (28.3) 
... ( ... ) 
... ( ... ) 
15.7 (52.6) t 
4.05 (28.3) t 



79.1 
14.0 
0.90 
24.7 
0.45 
2.25 
143. 
623. 
400. 
1242. 
1585. 
441. 

3978. 

2007. 



(231.)t 

(43.2) t 
(6.75)t 

(95.3) t 
(4.05)t 
(10.3)t 

(318.) 
(1174.) 
(1006.) 

(2670.) 

(3274.) 

(945.) 
(6335.)t 
(3518.) t 



Notes: total number of AGN inside detectable the SAFARI spectral range; f- 
outside the SAFARI spectral range. 



outside the SAFARI spectral range;!: excluding detections 
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TABLE 4 

Number, of Starburst galaxies detectable in a SAFARI survey of 0.5^ in IR lines as a function of redshift at So- {?,a) 

IN 1 TIR. IXTECRATIOX PER FoV. FOLLOWINC ORT'PnOXI ET AL. (2011) 



line/redshift 0<z<0.75 0.75<z<1.25 1.25<z<1.75 1.75<z<2.25 2.25<z<2.75 2.75<z<4 all redhifts # 

PAH(11.25//m) 515. (779.) f 245. (329.) f 270. (373.) f 71.1 (140.) 19.8 (48.2) 2.70 (10.4) 93.6 (199.) t 

[Nell] 12.81/im 326. (515.) f 212. (284.) f 128. (219.) 48.1 (102.) 11.7 (31.5) 1.35 (5.85) 189. (358.) j 

[Nelll] 15.55;im 62.1 (115.) f 52.2 (72.0) 29.2 (64.8) 13.0 (32.4) 2.25 (7.20) 0.45 (2.70) 97.1 (179.) j 

H2(17.03Mm) 26.6 (46.3)1 9.90 (24.8) 6.30 (11.2) 0.90 (4.05) ■•• (0.45) ••■ ( •■■ ) 17.1 (40.5) j 

SIII] 18.71/im 115. (199.) f 36.4 (72.0) 18.4 (43.6) 7.65 (20.7) 2.25 (7.20) • ■ ■ (0.45) 64.7 (144.) t 

O/y] 25.89/im 8.10 (15.3) 3.60 (9.90) 3.15 (11.2) 0.90 (1.80) ■•■ ( ■■■ ) ■■■ ( ■■■ ) 15.8(38.2) 

SIII] 33.48Mm 199. (326.) 180. (245.) 171. (270.) 71.1 (140.) 31.5 (69.8) 5.85 (18.4) 658. (1069.) 

Sill] 34.8lMm 256. (412.) 245. (329.) 270. (373.) 140. (240.) 69.8 (136.) 18.4 (46.4) 999. (1536.) 

Om] 51.8lAtm 92.7 (219.) 180.(245.) 171.(321.) 48.1 (140.) 4.05 (19.8) •■■ (0.45) 496.(945..) 

NIII] 57.32^111 153. (326.) 52.2 (122.) 18.4 (64.8) 0.45 (4.05) ■ • ■ ( • • ■ ) ■ ■ • ( ■ • • ) 224. (517.) 

OI] 63.18/im 637. (939.) 329. (459.) 321. (422.) 140. (240.) 31.5 (99.0) 5.85 (18.4) f 1458. (2159.) J 

OIII] 88.35;im 637. (939.) 212. (284.) 171. (321.) 71.1 (140.) f 19.8 (48.1) f 1-35 (10.4) f 1020. (1544.) j 

Notes: #: total number of AGN inside detectable the SAFARI spectral range; f: outside the SAFARI spectral range;J: excluding 
detections outside the SAFARI spectral range. 



TABLE 5 

Total number of galaxies detectable in a SAFARI survey of 0.5^ in IR lines as a function of redshift at 5a (Scr) in 1 hr. 

INTEGRATION PER FoV, FOLLOWING VALIANTE ET AL. (2009) 





line/redshift 


0<z<0.75 


0.75<z<1.25 


1.25<z<1.75 


1.75<z<2.25 


2.25<z<2.75 


2.75<z<4 


all redhifts # 


PAH(11.25Mm) 


3353 (5602) f 


2383. (4094.) f 


1074. (2743.) t 


322. (676.) 


91.3 (308.) 


4.95 (19.8) 


418. (1004.)t 


[Nell] 12.81/im 


1114 (1811) t 


251. (595.) t 


28.8 (136.) 


11.7 (49.1) 


1.8 (15.8) 


... ( ... ) 


42.3 (201. )t 


[Nelll] 15.55/im 


850 (1429) t 


94.0 (252.) 


51.8 (136.) 


26.1 (85.1) 


7.2 (29.2) 


■ ■ • (4.95) 


179. (507.)t 


H2(17.03/im) 


158 (331) t 


4.95 (17.1) 


1.35 (7.65) 


... ( ... ) 


... ( ... ) 


... ( ... ) 


6.3 (24.8)t 




SIII] 18.71;im 


158 (331) t 


4.95 (17.1) 


1.35 (7.65) 


... ( ... ) 


■■■ (1.80) 


... ( ... ) 


6.3 (26.5)t 




OIV] 25.89/im 


638 (1114) 


94.0 (252.) 


136. (327.) 


49.1 (211.) 


15.8 (51.8) 


■ ■ ■ (4.95) 


938. (1961.) 




5///] 33.48/im 


638 (1114) 


251. (879.) 


136. (327.) 


49.1 (211.) 


29.3 (142.) 


■ ■ ■ (4.95) 


1103. (2678.) 




5*//] 34.8lAtm 


1114 (1811) 


879. (1761.) 


505. (1074.) 


322. (676.) 


207. (452.) 


9.9 (61.2) 


3037. (5835.) 




OIII] 51.8lAtm 


2258 (4010) 


879. (1761.) 


505. (1512.) 


212. (676.) 


29.2 (142.) 


■ ■ ■ (4.95) 


3883. (8106.) 




NIII] 57.32^4111 


850 (1811) 


94.0 (392.) 


28.8 (136.) 


■■■ (11.7) 


... ( ... ) 


... ( ... ) 


973. (2351.) 




OI] 63.18^im 


3353 (4762) 


3149. (5211.) 


1512. (2743.) 


676. (1326.) 


207. (639.) 


19.8 (61.2) t 


8897. (14681.) J 




OIII] 88.35/um 


3353 (4762) 


1263. (2383.) 


505. (1512.) 


322. (955.) t 


142. (307.)t 


4.9 (36.0)t 


5121. (8657) t 



Notes: #: total number of AGN inside detectable the SAFARI spectral range; f: outside the SAFARI spectral range;$: excluding detections 
outside the SAFARI spectral range. 



TABLE 6 

Number of galaxies detectable in a CCAT survey of 0.5^ m IR lines as a function of redshift at 5a {3a) m 

1 HR. integration PER FoV, FOLLOWING FRANCESCHINI ET AL. (2010) 



line/redshift 


0<z<0.75 


0.75<z<1.25 


1.25<z<1.75 


1.75<z<2.25 


2.25<z<2.75 


2.75<z<4 


[OI] 63.18//m 


... ( ... ) 


... ( ... ) 


... ( ... ) 


0.27 (3.60) o 


• • • (0.27) t 


■■• (0.22)t 


[OIII] 88.35/im 


... ( ... ) 


22.2 (67.5) o 


1.80 (25.8) t 


2.07 (12.1) t 


86.0 (219.) # 


21.1 (67.0) $ 


[Nil] 121.90/im 


13.5 (39.0) t 


0.22 (0.85) t 


0.27 (1.80) t 


12.1 (59.0) # 


3.15 (27.0) $ 


0.22 (1.53) § 


[OI] 145.52/xm 


0.45 (1.35) t 


• ■ • (0.09) t 


0.86 (3.02) # 


■ ■ ■ (0.25) $ 


... ( ... ) § 


■ ■ ■ (0.72) 1 


[CII] 157.74//m 


235. (643.) t 


1786. (3951.) # 


1314. (1998.) $ 


517. (1116.) $ 


219. (603.) § 


350. (603.) 1 



Notes: o: at 200;Ltm band; f: at 230jum band; $: at 291;y.m band; #: at 350/Ltm band; $: at 450;Ltm band; §: at 620;y.m band; 
^:at 740/J.m band 



FIR spectroscopic cosmological surveys 



TABLE 7 

Number, of galaxies detectable in a CCAT survey of 0.5^ in IR lines as a function of redshift at 5<t (3(t) in 
1 iir. integration per fov. followixc grrrpioni et al. (2011) 



line/redshift 


0<z<0.75 


0.75<z<1.25 


1.25<z<1.75 


1.75<z<2.25 


2.25<z<2.75 


2.75<z<4 



O/] 63.18/im ...(...) ...(...) ...(...) 2.02 (13.7) 0.36 (1.53)1 0.22 (1.71) t 

OIII] 88.35;im • ■ • ( • • • ) 22.1 (53.5) o 5.49 (28.6) f 7.47 (25.3) t 108. (240.) # 39.7 (118.) $ 

Nil] 121.90/im 12.7 (32.5) f 0.81 (3.01) f 1-53 (5.49 t 25.3 (71.1) # 11.0 (37.9) § 1.71 (6.12) § 

OI] 145.52/im 0.81 (2.02) \ ■■■ (0.45) t 2.65 (9.85) # 1.00 (3.60) $ 0.13 (0.72) § 0.49 (3.01) 1 

CII] 157.74/im 172. (432.) f 1728. (3262.) # 1498. (2628.) $ 473. (1089.) $ 240. (534.) § 572. (882.) f 



Notes: o: at 200jLtm band; f: at 230;:im band; $: at 291;:im band; #: at 350jum band; $: at 450jum band; §: at 620iJ,ni band; 
t:at 740/:jm band 



TABLE 8 

Number of galaxies detectable in a CCAT survey of 0.5^ in IR lines as a function of redshift at Sct (So-) in 1 

IIR. intecration per FoV. followixc \'aliaxte et al. (2009") 



line/redshift 


0<z<0.75 


0.75<z<1.25 


1.25<z<1.75 


1.75<z<2.25 


2.25<z<2.75 


2.75<z<4 




OI] 63.18Aim 


... ( ... ) 


... ( ... ) 


... ( ... ) 


• ■ ■ (26.4) o 


■■ (1.94) t 


• • • (0.09) t 




OIII] 88.35//m 


... ( ... ) 


17.5 (54.9) o 


3.69 (28.7) t 


12.0 (49.5) t 


208. (450.) # 


36.1 (97.6) $ 




Nil] 121.90//m 


13.1 (40.5) t 


0.13 (1.00 t 


0.13 (3.69) t 


49.5 (136.) # 


29.2 (91.3) $ 


0.09 (4.81) § 




OI] 145.52/im 


0.94 (1.44) t 


•■■ (0.09) t 


1.35 (7.65) # 


■ • • (3.73) $ 


... ( ... )§ 


■ • • (1-66) 1 




CII] 157.74/im 


231. (639.) t 


2385. (5211.) # 


2061. (3586.) $ 


954. (2371.) $ 


450. (1215.) § 


531. (1075.) 1 



Notes: O: at 200;Ltm band; f: at band; J: at 291jum band; #: at 350;Ltm band; $: at 450/Ltm band; §: at band; ^:at 

740^m band 



